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Abstract. 

A normal form transformation is carried out on the operators of a complete 
set of commuting observables in a multidimensional, integrable quantum system, 
mapping them by unitary conjugation into functions of the harmonic oscillators 
in the various degrees of freedom. The transformation works at the level of the 
Weyl symbols of the operators, which are manipulated as formal power series in 
h by use of the Moyal star product. It is assumed that the Weyl symbol of one of 
the operators (the Hamiltonian) has a generic, stable fixed point in phase space. 
The normal form transformation takes place in a neighborhood of this fixed point. 
Once the normal form has been achieved, the Einstein-Brillouin-Keller or torus 
quantization rule follows easily, including higher order corrections in h. Crucial 
parts of the normal form transformation are not obvious generalizations of the 
one-dimensional case, nor is final quantization rule. The result raises some issues 
of differential geometry not found in the one-dimensional case. 



PACS numbers: 03.65.Sq, 02.40.Gh, 03.65.Vf, 02.20.Sv 

1. Introduction 

In a previous article (Cargo et al 2004), hereinafter referred to as I, we have developed 
a normal form transformation for a onc-dimensional quantum Hamiltonian whose Weyl 
symbol has a generic, stable fixed point in phase space. In particular, we showed that 
such a Hamiltonian can be transformed via unitary conjugation into a function of the 
harmonic oscillator in a neighborhood of the fixed point, in the sense of formal power 
series in h. One of the consequences is the Bohr-Sommerfeld rule for the eigenvalues 
of the Hamiltonian, extended to higher order in h. In this article we generalize 
this normal form transformation to multidimensional, integrable systems, in which 
a complete set of commuting observables, one of which has a generic, stable fixed 
point in phase space, is transformed via unitary conjugation into a set of functions 
of the harmonic oscillators in the various degrees of freedom. One consequence is 
the Einstein-Brillouin-KcUer (EBK) or torus quantization rule, including higher order 
corrections in h. Part of this calculation is a straightforward generalization of the 
methods of I, and part involves nontrivial new techniques. In the following a familiarity 
with I is assumed. References to equations, sections and appendices of I arc preceded 
by "I" (for example, Equation (1.17), Section 1.3.4, or Appendix LB). 

In our opinion the normal form transformation itself, which has many variations 
not yet explored (unstable fixed points, regions including separatrices, nongeneric 
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fixed points, etc.), and the diagrammatic and other techniques developed for carrying 
out the normal form transformation, are the most important results of I and of this 
paper. Nevertheless, the torus quantization rule and its higher order corrections are 
the results that are most easily appreciated and understood, so we have given them 
particular prominence in the following discussion. 

The usual torus quantization rule (Einstein, 1917; Brillouin, 1926; Keller, 1958; 
Percival, 1973; Berry and Tabor, 1976) is valid only for integrable systems, which at 
the level of the usual formula means that the classical Hamiltonian is integrable in a 
classical sense. If H is a classical, integrable Hamiltonian of N degrees of freedom with 
actions A = {A^, . . . , A'^) and angles (f> — {(p^, . . . , (p^), then the torus quantization 
rule states that the eigenvalues of the corresponding quantum Hamiltonian H are given 
approximately by £"„ = H{A), where H is expressed as a function of the actions, where 
the actions are evaluated at A^ = (nj + /A)h, where n — {ui, . . . , un) is a vector of 
integer quantum numbers, and where 7 = (71, . . . , jn) is the vector of Maslov indices 
(Maslov and Fedoriuk 1981, Littlejohn and Robbins 1987, Eckhardt and Wintgen 1991, 
de Gosson 1997, Foxman and Robbins 2004). See |Appendix A| for notation. The torus 
quantization rule is discussed more fully and precisely in the explanation of lO below. 

Examples of integrable systems include central force motion in any number of 
dimensions in ordinary particle mechanics, objects such as the 3nj-symbols in angular 
momentum theory (which are the eigenfunctions of complete sets of commuting 
observables), and integrable lattice and soliton systems. Integrable systems have 
also exerted a powerful influence on various areas of mathematics and mathematical 
physics, including symplectic geometry and representation theory (Weinstein 1977, 
Marsden and Ratiu 1999, Kirillov 1976). They are also the reference point for 
understanding partially integrable or nonintegrable systems. 

A normal form transformation of the kind developed in I has been studied by 
Helffer and Robert (1984) in the one-dimensional case. These authors use Fourier 
integral operators to carry out the normal form (rather than the Moyal star product, 
as in our work), and they prove the existence of the asymptotic H series for the 
eigenvalues. They do not, however, explicitly compute the higher order terms. In 
related work, Charbonnel (1986) has studied a multidimensional example. 

One question that arises immediately when contemplating higher order 
corrections to the torus quantization rule is what shall be the definition of an integrable 
system in quantum mechanics. A set of N operators need only commute through order 
h for their principal symbols to Poisson commute, and thereby to define a classically 
integrable system. The latter is all that is required for the usual (lowest order) torus 
quantization formula to be meaningful. When we extend torus quantization to higher 
order in h, is Poisson commutativity of the principal symbols sufficient, or must the 
operators commute to higher order in h? One might suppose that terms that are 
higher order in h could be handled by some kind of perturbation theory, and that 
classical integrability would suffice. 

In fact, it turns out that for torus quantization to be meaningful to all orders in 
h it is necessary to have a complete set of operators that commute to all orders in 
h. Without this condition, the normal form transformation that is the basis of our 
analysis does not exist, at least in the form of the power series in H that we posit for 
it. If this and other, reasonable (generic) conditions to be described below hold, then 
as we shall show the torus quantization rule can be extended to arbitrary order in h. 

We now present the generalized torus quantization rule, deferring the details 
of our assumptions to the main body of the paper. Let {H^, ■ ■ ■ , H'^) be a set of 
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N commuting, Hermitian operators, [H^,H''] — 0, j,k = 1, . . . ,N, oi which is 
"the" Hamiltonian (above referred to as simply H). Let the (Weyl) symbols of these 
operators be {H^ , . . . , H^) (that is, without the hat. See [Appendix A| ) These are 
assumed to have an S-expansion of the form, 

W = Hi + h^Hi + h^Hi + . . . , (1) 

that is, they are "slowly varying" or "admissible" (Voros 1977) with only even powers 
of h. The leading terms are the "principal symbols" of the operators . We 
assume that odd powers of h are absent from this expansion because their inclusion 
serves only to complicate all the formulas of this analysis without introducing any 
fundamentally new issues. In addition, most of the examples we can think of (such 
as {H,L^,Lz) in a 3-dimensional central force problem) have symbols that are even 
series in h, a property that evidently is related to the Hermiticity of the operators . 
The principal symbols Hq Poisson commute, {Hq, Hq} = 0, since all correction terms 
in the Moyal series (I. A. 6) for the symbol of a commutator are of higher order in h. 

The principal symbol Hq of the quantum Hamiltonian H^ will be treated as the 
classical Hamiltonian (above referred to simply as H). We assume that it possesses 
a generic extremum at some (fixed) point of phase space, so that classical motion in 
a neighborhood of this point is bounded and the fixed point is stable. Action-angle 
variables {(f>, A) are defined in this neighborhood, and the principal symbols can be 
expressed as functions of the actions, Hq = Hq{A). 

Then the eigenvalues of H^ are given by 



^io.kim {Af^ ^ A' ^ ^™) ) + Oih^] 

A={n+^/4)h 

where n and 7 are A^- vectors of integer quantum numbers and Maslov indices, as above. 
The ranges of the integers Uj depend on the contours on the tori used to defined the 
actions A^ , in a manner that will be made precise in subsection 12.31 Equation (0) is 
the torus quantization rule, carried through order h^. 

Equation jSJ) uses the following notation. (See also [Appendix A| ) Repeated 
indices are summed. The frequency matrix of the principal symbols and its action 
derivatives are denoted by 

_ dH^ _ dujjk _ d^ujjk , . 

These are all functions of the actions A. The angle bracket ( )^ at order in (|2J) 
indicates an average over the angles cj) (a torus average). The quantities inside the 
angle bracket depend on both A and 0, but of course only on A after the averaging. 
The diagrammatic notation is explained in Appendix LA, with obvious generalizations 
to the multidimensional case. The notation D'^ stands for the action-angle variables 
collectively, specifically, 

D'^ = i<l^\...,^^,A\...,A^) 

D,, ^ J^, D-- = {~A\ . . . ,-A^ ,cl>\. . . ,cl,^) 

(see [Appendix A)) . Finally, the entire expression is evaluated at the quantized values 
of the actions, A^ — {nj + jj)h. 
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Equation Q may be compared to its one-dimensional counterpart (I.l). The one- 
dimensional quantization rule at order fi.^ is expressed in terms of a single diagram, 
the second Moyal bracket {H, H}2 = H =5: H of H with itself. If H is expressed 
as a function of A (the action of the one-dimensional problem), this one diagram can 
be expanded by the chain rule yielding several diagrams, as in the multidimensional 
formula Equation (O, on the other hand, cannot in general be condensed into 
diagrams involving the symbols alone, but rather must be expressed in terms of 
diagrams involving the actions , and, in the case of the first diagram, also the 
angles (fP . It does reduce to (I.l) in the one-dimensional case, but it is not a simple or 
obvious generalization of the one-dimensional formula. Equation (0) will be discussed 
more fully in Sec.01 but for now we simply remark that of the three diagrams shown, 
the final two are in a sense "trivial" consequences of the formula for the symbol of 
a function of an operator, while the first one is the difficult one that requires new 
techniques for its evaluation. 

The normal form transformation of this paper takes place on three levels, of which 
the first two are purely classical. Normal forms of classical Hamiltonian systems have 
received considerable attention in recent years (for example, Eliasson 1990, Fomenko 
1991). The first level takes place in a small neighborhood of the fixed point of 
the classical Hamiltonian Hq, in which the functions Hq are represented by their 
quadratic approximations about the fixed point and a linear canonical transformation 
is used to achieve a classical normal form. (Specifically, it transforms these quadratic 
approximations into linear combinations of harmonic oscillators.) At the next level, 
this linear canonical transformation in a small neighborhood of the fixed point is 
extended to a nonlinear canonical transformation in a certain "region of interest," in 
which the level sets of the functions Hq are tori, achieving a nonlinear classical normal 
form in that region. (Specifically, the Hq are transformed into nonlinear functions of 
the harmonic oscillators). At the third level, the nonlinear classical normal form 
transformation is used as the skeleton for the construction of a quantum normal form 
transformation, valid in the sense of formal power series in H, constructed out of the 
Moyal star product. (Specifically, the symbols of the operators become functions 
of the symbols of the quantum harmonic oscillators in the region of interest). The 
quantum normal form then allows an easy derivation of the torus quantization rule. 

This paper is organized as follows. Section [21 contains an overview of the 
calculation, allowing readers who are familiar with I to understand the main ideas 
and how they differ from I. Section |31 contains a more detailed discussion of various 
steps in the presentation, including the rich subject of the classical normal form. 
Section 01 contains a discussion of various aspects of the result J^J, especially certain 
issues of differential geometry raised by the result. Section [3 contains conclusions. 
Finally, a series of appendices contain much of the real work of the calculation, as well 
as a detailed presentation of an example. 

2. An overview of the calculation 

In this section we present an overview of the calculation, assuming a familiarity with 
I, and we point out various similarities and differences with I. 

One difference concerns the notation for Weyl symbols and their expansions in 
powers of h. In I, the classical Hamiltonian was taken to be the symbol H of the 
quantum Hamiltonian H, that is, the whole symbol, possibly with an h dependence, 
whereas in this paper the classical Hamiltonian is taken to be the principal symbol 
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Hq of the quantum Hamiltonian H^, and similarly we use the principal symbols Hq of 
the "other Hamiltonians" together with iJg to define a classically integrable system. 
These principal symbols are independent of h. This approach is necessary in the 
multidimensional case because the full symbols of the operators do not Poisson 
commute in general, except when = 1 . A consequence is that in this paper notations 
for functions on phase space such as H^, K^, Gn, etc., are independent of h, whereas 
in I the analogous notations were allowed to have an h dependence. 

As in I, we begin by performing a canonical scaling of the operators Xj and pj 
to give them both units of action ' . This allows us to use a symmetrized form 
in Q below for the harmonic oscillator actions, and it means that all subsequent 
transformations preserve physical units. See [Appendix D] for an example. 

2.1. The two stages 

The main part of the calculation is to transform the set of operators {H^} via unitary 
conjugation into a new set of operators {M^} that are functions of the harmonic 
oscillators {/■' } in the various degrees of freedom. The harmonic oscillators and their 
symbols are given by 

I'^li^j+Pj)^ P{z)^\{x^,+p]) (5) 

(see [Appendix A| ). The or P are standard harmonic oscillators of unit frequency, 
which are their own actions. The new operators A'P are considered functions of the 
harmonic oscillators if their symbols A'P arc functions of the symbols in a 
neighborhood of the generic, stable fixed point of the classical Hamiltonian H^, in the 
sense of formal power series in h. We denote the assumed fixed point by zq- Here we 
rely on the theorem, the one-dimensional analog of which was discussed in I, that an 
operator is a function of the operators P if and only if its symbol is a function of the 
symbols (see [Appendix B| . 

As in I, the transformation from the operators to operators NP proceeds in 
two stages. In the first stage (the "preparatory transformation") we seek a unitary 
transformation [/ such that the operators K\ defined by 

= UHW\ (6) 
have principal symbols i^Q that are functions of the harmonic oscillator actions P . 
That is, the symbol of the operator is required to be a function of the harmonic 
oscillator actions {P} to lowest order in h. 

In the second stage, we seek a sequence of near-identity unitary operators C/2, IJa, 
Uq, etc., such that the operators AP, defined by 

= VWV\ (7) 

where V = . . . UQU4IJ2U as in (1.65), have symbols AP that are formally functions of 
the harmonic oscillator actions P to all orders in h. Each C/„ is given in terms of a 
Hermitian generator Gn by (1.47), and is responsible for transforming the symbols of 
the Hamiltonians into functions of the at order fi," (only even n are used). 
The two stages can be summarized schematically by 

{W} {k^} {AP}. (8) 

All this is very much as in I, but note that a single unitary operator U simultaneously 
transforms all the operators into the operators k^ , and a single set of operators 
{[/„} transforms all the operators k^ into operators AIK 



Multidimensional quantum normal forms, Moyal star product, and torus quantization^ 
Table 1. Notation for operators, symbols and functions depending on e. 
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2.2. The preparatory transformation 

The multidimensional method for finding U and the operators is similar to the 
one-dimensional case. First wc assume that U is given, and we work out expressions 
for the symbol of the transformed Hamiltonians . Then we impose conditions 
on these symbols (namely, that the principal symbols should be functions of the 
P). We assume that U can be imbedded in a smooth family Ue of unitary operators 
with boundary values at e = and e = 1 as shown in Table ^ The operators U 
of interest to us will satisfy this condition. Next we define the Hermitian operator 
Ge by (1.14), exactly as in I. The symbol of operator of interest to us will 
be constructed out of the classical mechanics associated with the integrable system 
defined by the principal symbols iJp . Since these have no h dependence, neither does 
Ge (unlike the case of I) . 

Next, as in (1. 11), we define 

Hi = Ulkm,, (9) 

so that Hi has values at e = 0, 1 indicated in Tabled Differentiating ©, we obtain 
a generalization of (1.13), 

^4p„A>,, (10, 

As in I, this equation is solved by transcribing it to symbols and expanding the solution 
HI in powers of h (see (1.17)), 

Hi = Hi, + h^Hl + h^Hl + . . . , (11) 

whereupon we obtain a set of coupled differential equations. 



DHL 
De 



= 0, (12a) 



-^ = -{G„i7i5J3, (12&) 
D77^ 1 1 

4e ^ ffl \ jC m\5, (12c) 

De 24^ 1920 °^ 

etc., exactly as in (1.16) except for the j superscripts. The bracket notation is explained 
in Appendix LA, and the notation D/ De is explained in Sec. 1.3.4. Similarly expanding 
the symbols , 

=Kl+ n^Ki + h'^Ki + . . . , (13) 

we find that the functions i/j^^ in l|ll|) are the solutions of the above differential 
equations subject to the initial conditions Hij^, — at e = 0. This implies H^, = H^ 
at e = 1, as shown in Tabled The boundary conditions at e = 1 differ from those in 
I, because of the manner in which the h dependence is handled. 
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The symbol is related to an e-dependent canonical transformation 

(symplectomorphism) : M^a^ _^ r2N ^^^^^ components Z^ : M^a' ^ (see 
Appendix LB for notation), exactly as explained in Sec. 1.3.3. As in I, when e = 1 
we drop the e subscript and write simply Z or Z^, as indicated in Tabled Then the 
solution of (|12a|) is 



H^>,^K^o°Z,, (14) 



or, at e = 1, 



Hi=KloZ, Ki=HioZ-'. (15) 

See (1.28) and (1.29). 

Now we choose Z so that is a function of the harmonic oscillator actions 
/■' . As in I, this is accomplished by using the action-angle variables (0, A) of the 
integrable system defined by the Poisson commuting functions and those defined 
by the set of harmonic oscillator actions {/^}, although the details are considerably 
more involved than in the one-dimensional case. The construction of the action-angle 
variables (0, A) is basically the standard one in classical mechanics, but one that for 
our purposes needs to be examined in some detail. This is done in subsections 13 . 11 - 13 . 31 
below. A consequence of this analysis is that over the region of interest, the functions 
Hq are invertible functions of the actions A' , 

Hi = fl o A, (16) 

where /q are functions : — > M, or, collectively, /o : M^. This is the obvious 

generalization of (1.35) (the 0-subscript means zeroth order in h, see subsection l2.3f) . 
Since the mapping /q is invertible in the region of interest, the frequency matrix ujjk 
of lO is nonsingular over this region. Another notation for the frequency matrix is 

ojjk = (17) 

(See Appendix LB for notation.) 

As for the harmonic oscillators these are their own actions. We define 
9 = {9^ , . . . ,0^) as the corresponding angles, so that {9^,P) is just the usual set 
of one-dimensional harmonic oscillator action-angle variables in the {xj,pj) plane. 
Then the canonical transformation Z satisfies 

A^ = P o Z, 

(18) 

(jyi ^9^ oZ, 

where for simplicity we have set the matrix v discussed in subsection l3.2l to the N x N 
identity. There is no loss of generality in this, since the conventional basis contours on 
the i?o-tori could be chosen to be the pull-backs of the standard contours on the I- 
tori. At the end of the calculation it is easy to carry out a linear transformation on the 
actions to restore an arbitrary choice of basis contours on the ifo-tori. Equation (|18|l 
is the obvious generalization of (L37) in the one-dimensional case. 
Now we combine the first half of (|18l) with (jKil) to obtain 

HloZ-'=Kl = floI, (19) 

where / denotes the map : with components . This shows that the Kq 

are functions of the P , as desired. 
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Having found Z , we imbed it in a family as in I, with values at e = 0, 1 as 
shown in Tabled and then we define intermediate (for any e) action-angle variables 

by 

Al ^ P o Z„ 

(20) 

cj^i ^e^o z,, 

as in (1.43). We also introduce collective notation for the various action-angle variables 
(see also (0}, 

Di^^{cj^l...,cj^^,Al...,A^), (21) 
i?^ = (^\...,0^, /^), (22) 
as indicated in Tabled Then combining H2()|l with (|14|l . we have 

H'oe = fi°A,, (23) 
as in (1.44). Then we define Ge as the solution of (L19), as explained in Appendix I.E. 
The inverse Weyl transform of is the operator Ge, which defines 11^ as the solution 
of (1.45), with value t/ at e = 1. 

This completes the solution of H12a|l and the construction of other objects needed 
for the 0-th order part of the preparatory transformation. The higher order terms are 
obtained as in I, for example, we integrate (jl2&|) to obtain 

(24) 



Ki = Hi oZ-^-^ J^' de {Ge, H^Js o Z;\ 



with an extra term compared to (1.32) because of the different boundary conditions 
at e = 1 (see Table [ij. 

2.3. The second stage 

The second stage transformation proceeds very much as in I, as indicated by (0 and 
(1.47). In particular, writing 

AP = + ffMi + h'^Ml + ... (25) 
for the expansion of the symbol of IS/P , we find the hierarchy of equations, 

A^=K^, (26a) 
Aq^Ki+{G2,K'o}, (266) 
Mi ^Ki + {G2,m ~ ^{G2,Klh 

+ i{G2,{G2,A'^}} + {G4,Xi5}, (26c) 

etc., just as in (1.51) except for the j superscripts. As in I, we wish to choose the 
generators G„ so that the are functions only of / (i.e., they are independent of 9). 
At zeroth order we have this already, 

M^=K'a^f^oI = H'^oZ-\ (27) 

as in (1.52). At second order we obtain by projecting out the averaged part of 
ll26^ . 

Aq^Ki, (28) 
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as in (1.54), where the overbar indicates an average over the angles 9 (a torus average). 
Projecting out the oscillatory part, we obtain a differential equation for G2, 

^ = -i^-%^Ki, (29) 

as in (1.55), where the tilde indicates the oscillatory part of a function and where we 
employ the same abuse of notation mentioned beneath (1.52). 

Unlike the one-dimensional case, H29|l is not guaranteed a solution unless an 
integrability condition is satisfied. This question is discussed in [Appendix E| where it 
is shown that the generators Gn exist to all (even) orders. The proof relies critically 
on the fact that the original operators commute to all orders in h. Without this 
condition, the normal form transformation does not exist as a power series in h. 

From this point it is straightforward to proceed as in I to obtain a generalization 
of (1.56), an expression for M| containing an e-integral of a third order Moyal bracket. 
This integral can be done, but it ends up reproducing some of the work involved in 
calculating the symbol of a function of an operator, a problem that is best handled by 
the methods of Gracia-Saz (2005). To take advantage of those methods we proceed 
somewhat differently. 

Granting that the generators Ge, G2, G4, etc. exist, we use the unitary 
transformations U2, etc. to transform the operators P into new action operators 
we denote by , reversing the sequence of steps used to transform into JVP as 
shown in {Tj) and (jSJ. That is, we define 

=^ ft ji V, (30) 

as in subsection 1.5.2. We call the "action operators," but we do not denote them 
since by our conventions the latter would be the Weyl transforms of the classical 
actions . The operators A^ are unitarily equivalent to the harmonic oscillator 
actions /-' only through lowest order in h, while the operators B^ are equivalent to 
them to all orders in h. Consequently, the operators B^ have symbols which are equal 
to the classical actions A^ at lowest order, that is, if we write 

B^ ^Bl + h^Bi + . . . , (31) 

then _Bq = AK As for the first correction term, we find 

(^ + ^^'rfe{G„^n3oZri)oZ, (32) 

where we abuse notation in the first term (see Appendix LB) by treating G2 both 
as a function of {9, 1) and as a map : R. The derivation of (|32|l is given in 

subsection l3.4l 

The transformations we have constructed guarantee that the operators are 
functions of the operators P , 

£p^f{i), (33) 

for some functions : R^ ^ K, or, collectively, / : R^ ^ R^. Equation ^ is 
the normal form for the original operators = A'PV , and it serves to define the 
functions . Conjugating by V and using {T)) and l|Sn)l gives 

W = V^f{i)V = fiV^iV) = fiB), (34) 

an equation which determines the eigenvalues of the once and the eigenvalues 
of the B^ are known. But by (|30|l the latter are the same as the eigenvalues of the P , 
namely, {uj + l/2)fi, where rij G N = {0, 1, . . .}. 
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As for the functions , to find them we take the symbols of both sides of H34|l , 
using the methods of Gracia-Saz (2005) for the symbol of f^{B). This gives 



-0(n*),(35) 



where and its derivatives are evaluated at B (more properly, they stand for 
o B, etc.) Then we expand H and B as in ^ and H31|l and similarly write 
= /q + + . . ., expanding H35|l and collecting terms by powers of h. At lowest 

order we have Hg = /o(^), showing that /g is the classical action-energy map (it is 

the same function /q introduced in (|16|l 'l. At second order we find 

m = fi+U.k B'2~^Um ^ ^0 +^/o.fc^™ i^' ^A^^ A^) , (36) 

where now and their derivatives are evaluated at A. We consider everything in this 
result to be functions of ((/>, A) and we take the average of both sides to obtain 

-YA^UimiA^'^A^^A^A , (37) 

/ 

where the angle brackets indicate an average over as in 

Equation (|37|) gives in terms of functions that have been determined, but it 
would be more satisfactory if the e-integral contained in the term B2 were done (see 
H32|l '). Notice that for the computation of /2 (and hence the quantization rule) we 
only need the (/)-average of Bt^- This average kills the first term in 1)32(1 . since 

See (|18|l and Appendix I.D. As for the e-integral (the second term in (|32(l '). it can be 
evaluated as follows: 



de (G, 3 Ai) oZ-^oZ^- A^ i^D^ + (39) 

where A" is a function that is not needed for the evaluation of fi^ or the quantization 
condition, since its ^-derivative vanishes upon 0-averaging. This result is derived in 
[Appendix F| Thus, we have 

{b',)^^1{d^^A^^D,)^. (40) 

Using now (|^ in P7| and writing ujjk — fo k ° A (see © and (|T71l ). we obtain 
the second order term of apart from the special choice of contours P , which 
forces the ranges of the quantum numbers nj and the Maslov indices to take on their 
values for the harmonic oscillator {nj E N and 7^ = 2). But the entire expression is 
covariant under a GL{N, Z) transformation of actions and angles, so if we set A' — vA, 
(p' = v~^(l) where v £ GL{N,Z), then we have n'j = J^k'^jkrik and 7^ — ^J^k'^jk- 
Dropping the primes, we then have the general formula Further discussion of this 
result is presented in Sec. ^ 
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3. Details of the calculation 

In this section we present a series of digressions from the overview of Sec.|21that supply 
details of the calculation. 

3.1. Classical normal form in a small neighborhood of zq 

In this subsection we develop the classical normal form transformation Z in a small 
neighborhood of the fixed point zq of the Hamiltonian Hq, in which the functions 
Hq are well approximated by their quadratic expansions about zq. (Properly, this 
discussion should be framed in terms of tangent spaces and germs of smooth functions, 
but we shall proceed somewhat intuitively.) At this level of description, Z is a linear 
canonical transformation (the linearization of the fully nonlinear Z defined over the 
region of interest). To begin we list the assumptions we make about the behavior of 
the principal symbols Hq in a neighborhood of the fixed point zq. These amount to 
taking the most generic case of a stable fixed point of an integrable system. 

We assume that all symbols (that is, all the terms of their h expansions) 
are smooth in a neighborhood of zq. This excludes the case of singular potentials, 
such as the Coulomb potential. We define as the 2N x 2N real, symmetric Hessian 
matrix of the principal symbol Hq evaluated at zq, 

- ^('.)^ (41) 

We assume that (the Hessian of the Hamiltonian) is positive or negative definite, 
corresponding to stable motion about the fixed point. For convenience, if is 
negative definite, we agree to replace by —H^ at the beginning of the calculation, 
so that in the following we can assume that is positive definite. The for 
j > ! need not be definite (positive or negative). We assume that the set of matrices 
{Q'^tj = l)---)-^} is linearly independent. For simplicity, we assume that all the 
principal symbols vanish at zq (if not, we just subtract constants from the operators 
to make this so). 

These assumptions imply that the first derivatives of the Hq at the fixed point, 
denoted by 

V = j;^(-o), (42) 

all vanish. For j = !, this follows from the assumption that zq is a fixed point of 
Hq. For j > !, we take the derivative of the vanishing Poisson bracket {Hq, Hq} and 
evaluate it at zq: 

d_ 

dzf" 



The second term on the right hand side vanishes since V^^ = 0, which implies that 
— for j > !, since both Q^^, and J^" are nonsingular matrices (in the former 

case, because is positive definite). Thus, the Hamiltonians Hq all have a common 
fixed point at zq, and their Taylor series expansion through quadratic order is 

H'oi^.) = l-Q^^,{z-zonz-zor. (44) 
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Now computing the second derivatives of the vanishing Poisson bracket {H^ , Hq} 
and evaluating at zq, we obtain the matrix equations, 

[JQ' , JQ''] ^ 0, (45) 

or 

Q^JQ'' = Q'^JQ^ = -{Q^JQ^f = antisymmetric, (46) 

where T indicates the transpose. 

Regarding the set of matrices {Q-*}, we have the following theorem. Write the 
harmonic oscillator actions in the form, 

P{z) = \{x^^+p]) = \z'^^^z, (47) 

which defines A-' is the 2N x 2N symmetric matrix that contains zeros everywhere 
except for ones in the (j, j) and [N+j, N+j) positions. Then there exists a symplectic 
matrix S such that 

S^Q^S = ajkA^ (48) 

where ajk is a real, N x N, nonsingular matrix. Moreover, the numbers aik (the 
first row of the matrix) are all positive. The proof of this theorem is discussed in 
[Appendix C| The proof relies on the fact that is positive definite. 

Now we define a canonical transformation Z in terms of the symplectic matrix S 
produced by this theorem, 

Z{z)^S-\z-zo). (49) 

This allows us to combine (|44ll and (|48|l to give 

H'o = ajk o Z. (50) 

Thus, under the transformation Z (or Z^^), the Hamiltonians Hq are mapped into 
functions of the harmonic oscillators P (in fact, linear functions). Therefore we 
identify this Z with the transformation Z used in H15|l . and we have the desired 
normal form transformation in a small neighborhood of the fixed point, 

A'i5=a,fe/^ (5f) 
3. 2. Tori, contours and actions in a small neighborhood of zq 

The level sets of the harmonic oscillator actions P (say, P = Cj, for Cj > 0) are n-tori, 
where n is the number of the Cj that do not vanish. Thus, < n < A^, and the subset 
of that contains tori of dimension < n is the union of a set of N\/n\{N — n)\ 
planes of dimension 2n passing through the origin. On the other hand, 15Uf) shows 
that Z maps the level sets of the functions Hq near zq onto those of the harmonic 
oscillators near the origin. Thus, the level sets of the near zq are also tori of 
dimension ranging from to that fit together exactly like those of the harmonic 
oscillators near the origin. We shall call these two sets of tori the "iJo-tori" and 
the "/-tori." 

To visualize how the tori of various dimension fit together, we consider the two 
maps, 

iJ„:R2^->R^:z^(i/i(z),...,<(z)) 

/:M2A^ ^R^:zk^ (/1(z),...,/^(z)) ^ ^ 



Multidimensional quantum normal forms, Moyal star product, and torus quantization!?! 



where can be thought of as "constant-of-motion space." These are examples of 
the energy-momentum map discussed by Cushman and Bates (1997). We call these 
the "i?ci-map" and the "/-map," respectively. These maps and others associated with 
them are illustrated in Fig.^ An entire level set of the functions iJg or is mapped 
by these two maps into a single point of constant-of-motion space. In the following we 
think of the domain of the i/o-map as a small region U around zo (an open ball) . The 
domain of the /-map may be taken to be all of M^^, but we shall often be interested in 
the small region V — Z{U) around the origin (so V is also an open ball). The images 
of these maps are the sets of values the constants of motion, Hq or P , can take on. 

The image of under the /-map is the region Cj > 0, for j — 1, . . . , N, for 
example, it is a single octant of M'^ when N = 3. Thus, the image of V under the 
/-map looks like a corner of an iV-cube that has been sliced off in such a manner that 
the cut does not pass through the vertex itself, as illustrated by the shaded area in 
the upper right of Fig.Q] The sliced corner has flat faces of dimensionality n, ranging 
from to TV, where the n = face is the vertex (the origin of K^), and the n = N 
face is the entire sliced corner. In the figure, the 1-faces are heavy lines, coinciding 
with the positive axes. The inverse image of a point in the interior of an n-face of the 
sliced corner is a torus in of dimensionality n. 

By throwing out all the faces with n < N, we can create an open, contractible 
region of K^, whose inverse image under the /-map is the set of invariant iV-tori of 
the harmonic oscillators. This set is topologically the Cartesian product of the just 
mentioned region and an iV-torus. Thus, the subset of phase space containing 
the iV-tori is connected and has the form of a trivial fiber bundle in which the fiber is 
an A''-torus. 

Now let Hq ~ for j 1, . . . , iV be a level set of the //q near zq, which is 
mapped by Z onto a level set /-' — Cj near the origin. Then by (|50() we have 

= Ojk Ck, (53) 
which means that the image of U under the //g-map is just a linear transformation 
(specified by ajk) applied to the image of V under the /-map (the sliced corner of 
the iV-cube). In particular, all the level sets of the functions Hq that are A^-tori 
correspond to interior points of the sliced corner of the iV-cube, while tori of lower 
dimension reside on the surface. Thus, the image of U under the //o-map has a 
"corner." See, for example, Fig. 3.3 in Chapter IV in Cushman and Bates (1997) (in 
that figure, the corner does not occur at the origin of M^, because the constants of 
motion do not take on the value at the fixed point). 

Now we consider basis contours on the various tori. A set of closed contours 
on a torus constitutes a basis if any closed contour can be represented as a linear 
combination of the basis contours with integer coefficients, with the understanding 
that two closed contours are considered equal if they can be continuously deformed 
into one another. That is, the contours form a basis of the first homology group H^{Z) 
of the torus, with homologous contours considered equal. On a specific Af-torus, the 
choices of basis contours are related by GL{N,'Z) transformations (see [Appendix A| |. 

As for the /-tori, the functions /■' (z) are themselves actions of these tori, that is, 
there exist basis contours {C^} on the /-tori such that 

I'i^) = / Pkdxk. (54) 

Specifically, we shall take the contour to be the circle in the Xj-pj plane whose 
radius is determined by the contour value of P , with the other Xk and pk fixed. This 
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gives a smooth assignment of contours on the /-tori over all of R2A'^ and hence 
over the region V. 

As for basis contours on the TJo-tori, we could take these to be Z^^ o , where 
we think of the as maps from the circle onto phase space, : M^^, but in 

practice other basis contours on the i/o-tori, call them F-', may be more convenient 
(for example, see [Appendix D| |. On each 7/o-torus, the contours must be related 
to ^ o by an element ly e GL{N, Z), so that 

ZoP = jyj-feC'=. (55) 

By requiring that i> be the same for all _ffo-tori, the contours constitute a smooth 
assignment of basis contours over the 7?o"toi'i in the small region U. 
Once the are given, the actions are defined for z G U 

A'{z) = -^j> Pkdxk, (56) 

and they are smooth functions over U. Then since the symplectic 2-form d{pkdxk) 
is invariant under the canonical transformation Z, we can map the integral H5f)(l from 
the 7/o-torus to the /-torus to obtain 

= i^jk I'' o Z, (57) 

which is the first half of H18(l . here established over small neighborhood U around zq. 

As for the angles, we first define the harmonic oscillator angles 6^ as the clockwise 
angles in the Xj-pj plane, vanishing on the positive Pj-axis (a convenient definition). 
Then the second half of H18|) can be taken as the definition of the angles 0^ of the 
//o-tori. 

Finally we combine H5()|l and (|57|l to obtain, 

//i^' - (al/-l)JfeA^ (58) 

which gives the form of the functions H^{A) in neighborhood U. Differentiating this 
and using Q, we obtain 

ujjk{0) = {au~^)jk, (59) 

where the frequency matrix is considered a function of the actions A, and is here 
evaluated &t A = Q (corresponding to the fixed point). This matrix is nonsingular, 
since both Ojk and z/j^ are nonsingular. 

This completes the classical normal form transformation in a small neighborhood 
of the fixed point zq. [Appendix D| contains a physical example of this transformation, 
in which all the constructions of this subsection are carried out explicitly. 



3. 3. Classical normal form in the region of interest 

In this subsection we discuss the definition of the "region of interest" , denoted R, an 
enlargement of the small neighborhood U of the fixed point zq which is the domain of 
the (generally nonlinear) canonical transformation Z that brings about the classical 
normal form (|15|l and satisfies (|18|l . The situation is illustrated in Fig. [3 which is a 
modified version of Fig.^ The region R and the fixed point zq are shown in the lower 
left corner of this figure. The transformation Z maps R onto the region W = Z{R), 
upon which the harmonic oscillator actions / are defined. The mappings Hq and / 
are the same in Fig. (21 as in Fig. except for the expansion of the domains. Note 
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that in subsection l3.1l the notation Z stood for the hnearization of the Z that we are 
discussing in this subsection. For example, H49|) must now be understood as 

— (zo) = (5-i)%, (60) 

with Z{zq) ~ 0. 

A simple approach to understanding the transformation Z is just to say that it 
maps the action-angle variables of the Hamiltonians into those of the harmonic 
oscillators , as in H18|l . This approach requires the action-angle variables for Hq 
to be defined. Then Z is canonical because it is the composition of two canonical 
transformations. Moreover it is easy to see that Z is smooth, except possibly where 
the A'^-tori degenerate into tori of lower dimension and some of the angles become 
undefined. In I worries about the smoothness of Z at the fixed point led us to take 
another approach, in which the Morse lemma was used first to construct Z, and then Z 
was used to define the angles (j) (see Appendix I.C). Nevertheless, the simpler approach 
is still useful, in that it shows that the domain of Z in the one-dimensional problem 
can be extended out from the fixed point to the first separatrix, where the action-angle 
variables become undefined. 

In the multidimensional case, the construction of Z can again be based on 
symplectic generalizations of the Morse lemma. This has been done by Eliasson (1990), 
in more generality than we need here. In the following we shall simply assume the 
existence of a smooth, symplectic map Z : R ^ such that Z{zq) = and such 
that ifo-tori are mapped onto /-tori. As in I, we shall take R to be the largest region 
containing Zq over which the entire construction works; R is required to be a union of 
invariant i?o-tori. In the multidimensional case, however, obstacles to the extension 
of the domain R include not only separatrices but also pinched tori and other kinds 
of nongeneric level sets of the functions Hq. The following remarks are useful for 
visualizing the kinds of such level sets that can arise in the multidimensional case. 

First we comment on the Liouville-Arnold theorem, which says that if the 
differentials dHQ are linearly independent on a compact, connected component of 
a level set of the Poisson commuting functions Hq , then that level set is a Lagrangian 
A^-torus. By the assumptions we have made, the levels sets of the Hq are certainly 
compact in a sufficiently small neighborhood of zq, because the level sets of Hq alone 
are compact there (they are the spheres S'^^~^). But the differentials dHg are not 
necessarily linearly independent on those level sets, even in a small neighborhood of 
Zq, as shown by the invariant n-tori for n < N, which exist arbitrarily close to zq 
in the manner described in subsection 13.21 For example, at the fixed point itself (a 
0-torus), the differentials dHQ all vanish. 

We consider the rank of the set of differentials {dHQ} to be the rank of the 
N X 2N matrix Afj^ — OHq/Oz^. This rank is constant along any Hamiltonian flow 
generated by any of the Hq , as are the Hamiltonian functions Hq themselves. The 
commuting Hamiltonian flows associated with the Hq deflne an R^-action on phase 
space, assuming that the advance maps exist for all time. Thus a level set of the 
Hq consists of a union of orbits of the R^-action, upon each of which the rank of 
the differentials dHQ is constant. Each of these orbits is diffeomorphic to R'' x {T^Y, 
where r + s = n (a generalized cylinder), where is the 1-torus (a circle), and where 
the dimension of the orbit n equals the rank of the differentials dHQ . 

The differentials dHg can fail to be linearly independent on a level set in different 
ways. The n-tori for n < N that exist even near the flxed point are composed of 
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a single orbit of the M^-action, upon which the rank of the differentials is n. The 
separatrix in a double well oscillator in one dimension consists of three orbits of 
the M^-action, the two orbits that asymptote to the fixed point ds, t ^ ±c», each 
diffeomorphic to R, and the fixed point itself, a 0-torus. Generalizations of this 
behavior occur in higher dimensions. Another example is offered by the pinched 
torus of the spherical pendulum, which exhibits the phenomenon of monodromy 
(Duistermaat 1980, Cushman and Duistermaat 1988, Cushman and Bates 1997, Alber 
et al 1998, Child 1998, Cushman et al 2004). In this case the pinched torus minus the 
pinch point is a 2-dimensional surface diffeomorphic to M x S'^ , while the pinch point 
itself is a 0-torus. 

Another issue concerns the definitions of the actions A' as smooth functions 
: R M.^ , which requires a smooth assignment of basis contours over the i?o-tori 
in R. Consider a region of phase space in which the differentials d7?Q are linearly 
independent and the level sets of the Hq are A^-tori, and then consider the image 
of that region under the HQ-msq). This image can be taken to be the base space 
of a discrete, principal fiber bundle, in which the fiber over a point is the set of all 
basis contours on the corresponding A^-torus. The structure group of this bundle is 
GL{N,'Z), which relates the sets of basis contours. A smooth assignment of basis 
contours as a field over the base space is equivalent to a global section of this bundle, 
which exists if and only if the bundle is trivial. But such bundles are nontrivial 
precisely where monodromy occurs (this is the meaning of monodromy). Notice that 
if we have a smooth Z that maps _ffo-tori in some region onto /-tori, we can pull 
back the basis contours of the /-tori, which can be smoothly assigned, to create a 
smooth assignment of basis contours on the i/o-tori. Thus, monodromy cannot occur 
in the domain of any such Z, such as the region R considered here. In the example of 
monodromy given in Fig. 3.3, Chapter IV of Cushman and Bates (1997), one would 
have to introduce a branch cut terminating on the pinch point in order to create a 
region R free of monodromy. There is no unique way to do this, and hence no unique, 
maximal set R upon which Z is defined. 

Finally, we impose the requirement that R be contractible, since this guarantees 
the existence of the generator (see Appendix I.E). Note that this condition by itself 
guarantees that any bundle over R is trivial. 

To proceed we assume Z is defined on some suitable set /?, and that 
its linearization at zq agrees with the transformation previously worked out in 
subsection l3.1l A smooth assignment of basis contours F-' on the //o-tori in R is made 
according to ()55|l . where ly is the same matrix constructed in the small neighborhood 
of Zq (it is extended to all of R) . Then (|57|l holds in all of R, and maps A and in 
Fig. 13 are defined and satisfy A — o I o Z . One nonobvious conclusion from this 
is that the boundaries of the image of the A-map are flat (straight lines in Fig. ^ . 
Another conclusion is that the actions A^ uniquely label the //o-tori, that is, no two 
i/o-tori in R have the same actions. This follows since the actions P uniquely label 
the /-tori. 

This means that the map /o in Fig. [21 is defined, for if we take a point in the 
image of the A-map, its inverse image is an invariant //o-torus, upon which the 
values of the Hq are defined. Thus, Hq = fo ° A as in l(TB|) . Moreover, the map 
/o is of maximal rank (N) and therefore invertible in the region containing A^-tori, 
where the differentials dHg are linearly independent. This means that the frequency 
matrix ujjk{A) is invertible in this region. The frequency matrix can be defined on the 
boundaries of the image of the A-map by a limiting process; this limit is nonsingular 
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at the fixed point (see H59|l 'l. We shall assume that the limit at other points on the 
boundary is also nonsingular (this excludes the possibility that the kind of nongeneric 
behavior seen in the quartic oscillator in one dimension might develop as we move 
along the space of tori of dimension n < N). Thus, the image of the HQ-map is a 
smooth deformation of the truncated corner of the A^-cube, as illustrated in Fig. |2] 
(with corners and curved faces). These features are clear in the example of Fig. 3.3, 
Chapter IV of Cushman and Bates (1997). Finally, we define the map Kq as Hq o Z~^, 
as in H15|l : this implies Kq ^ fa o v o I , as in p9|l (where we have replaced v by the 
identity) . 

5.^. Symbols of the action operators 

The derivation of the symbols B^ of the action operators B^ , defined by (|3()|l . is 
a straightforward application of the method used to derive the symbols A'P from 
the symbols , except that it is run backwards. For convenience we define an 
intermediate set of operators by B^ = W U, that is, 

= U^UI...P ...U4U2. (61) 

Then by expanding the Un in terms of their generators and transcribing to symbols, 
we obtain 



= P +h'^{G2,P} + n'^ 



i{G2,{G2,P}} + {G4,P} 



0(^6), (62) 



where third and higher order Moyal brackets involving P are dropped since is a 
quadratic polynomial. 

Next we define B{ — , take e-derivatives to obtain a differential equation 

for the B^ , and transcribe them to symbols, exactly as in the steps leading up to (|12afl - 
p2c|l . The ^.-expansion of the symbols B( yield terms B-j^^ with boundary values at 
e = 0, 1 shown in Tabled Note that Cq = The differential equations are easily 
solved, yielding Bq — and (|32|) through second order. 



4. Discussion of the torus quantization rule 

In this section we discuss various aspects of the torus quantization rule ^ , the normal 
form transformation, and the calculation used to derive them. 



Doing the e-integrals 

In both I and this paper, the e-integrals that arose in the computation of the normal 
form function / all turned out to be doable, at least through second order in h, 
although the methods employed were ad hoc and no proof was given that the same 
could be done at higher orders. See Appendix I.G and [Appendix F] In principle it 
would be possible to carry out the normal form transformation using the integral 
representation of the Moyal star product, so that the conjugation © could be 
evaluated by the stationary phase approximation. This approach would avoid the 
necessity of imbedding U ot Z in a one-parameter family of transformations, and 
there would be no e-dependent paths through various group manifolds to integrate 
over. (But one would have to work with the rapidly oscillating symbol U of U, 
something we have avoided in the approach of this paper.) The existence of this 
alternative approach suggests that all e-integrals that arise in the present approach 
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that are needed for the normal form are doable, that is, the integrals can be expressed 
in terms of the transformations at the end points of the path (that is, evaluated at 
e = 0, 1), and do not depend on the path itself. An e-integral that is independent of 
path must mean that in some sense we are integrating an exact 1-form on the infinite- 
dimensional group manifold. See Appendix I.H, in which the e-integral occurring 
in the expression for K2 in the one-dimensional case was shown to change by an 
exact 0-derivative upon a variation of the path, that is, exactly what was needed to 
guarantee that the normal form be invariant. The calculation of the variation of the 
e-integral in that case relied on the star version of the Jacobi identity. It is likely 
that a general understanding of the path dependence of the e-integrals in the general 
case (multidimensional, arbitrary order in h) would result by following the suggestions 
raised by this calculation. 

4-2. The quantization rule involves the angles 

A striking aspect of the quantization rule ^ is the presence of the angles angles 0'^ 
in the first diagram of the second order term, which we see if we write D^^ = (0, A). 
This diagram is the difficult one derived in [Appendix F| The appearance of angles in 
the quantization condition requires some discussion. 

First, we note that in the one-dimensional problem the angles can be eliminated 
from the quantization condition. This is shown to second order in h by the actual 
result in I, and to fourth order by Gracia-Saz (2005), who used Colin de Verdiere's 
(2004) method for computing the eigenvalue series combined with his own theory of 
the symbol of a function of an operator. In fact, it is implicit in Gracia-Saz's work 
that the angles can be eliminated to all orders in h in the one-dimensional problem. 

Next, we note that the angles are not unique, since their origin on each torus 
(the point where the angles are taken to vanish) can be chosen with a great deal of 
freedom. That is, given one set of action-angle variables {(j),A), it is always possible 
to subject them to a canonical transformation of the form. 



where F{A) is some function of the actions. Geometrically, the surface = 
is required to be a Lagrangian manifold. The quantization condition must be 
independent of such a shift in the origin of the angles. 

To show that it is, we write the 0-dependent diagram in Q as 



Under the replacement (j)^ 1-^ cj)'' + F^k{A), the right hand side goes into itself plus 
the sum of three terms, of which two cancel by simple symmetry considerations. The 
remaining term is 




(63) 



A'J = A^ 



^A^^Df, = - 0fc ^ ^3 =S A*^ + A^" Ai =S 0' 



,fe 



(64) 



A' 



F^klm A'' A^ 



(65) 




j^i ^ A™)] , (66) 
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which vanishes upon the (/>-averaging in |(2Jl. Therefore, although the diagram (|54|l 
itself is not invariant under the shift of origin l|63|) . its angle average is, and hence so 
is the quantization condition 

A third aspect of the appearance of angles in the torus quantization rule (|2Jl is 
that in practice angles are much harder to compute than actions, either numerically 
or as explicit functions of cc's and p's, and what is worse is computing the derivatives 
of the angles with respect to z^. This makes a numerical test of the quantization rule 
^ difficult. Since the quantization condition itself does not depend on the origin of 
the angles, one wonders whether there is a method of computing the second order 
term that does not rely on this origin, either. We do not know the answer, but the 
following are some relevant considerations. 

First, not all diagrams with angles cause difficulties. For example, the angle 
average of a diagram of the form (^'^ ^ (X) , where X is any diagram, can be written 
in the form. 




in which the operator d/dA^ has been taken out of the average. This operator, when 
applied to a function of both A and cf), is not invariant under the transformation Ht)3|) . 
but it is so when applied to a function of A only. In fact, such action derivatives, 
outside an averaging operator, occur in the one-dimensional quantization rule. 

One might hope to transform the difficult diagram into another one that is 
independent of the origin of the angles, plus terms that average to zero. We doubt 
that this can be done, but the following identity is a interesting attempt along these 
lines: 

D^^A'^^Df, = Dt'^iA^^D^) + / \ 



= -M^^'-^'^ -w^""'-^'^ ^''^ 

where to prove the first step we expand the first term after the equals sign by the chain 
rule and use the expansion of the vanishing diagram (D^ — Z?^) , and where 

the second step follows by using Q), expanding Poisson brackets in the coordinates, 
and using the symmetries of the triangle diagram (sec (|H.2p ) and (I. A. 3). The first 
term on the right hand side of H68|l is of the type indicated in l|67|) (hence, independent 
of origin of the angles after averaging) and the second term vanishes upon averaging. 
Unfortunately, the third term still depends on the origin of the angles, except in one 
dimension (see subsection I4.4|l . 

4-3. Changing the independent variables 

By introducing a symplectic connection (Bayen et al (1978), Fedosov (1994)) the 
derivatives with respect to z'' that appear in the definition of a diagram (see 
Appendix LA) can be replaced by covariant derivatives that can be evaluated in any 
other coordinates, of which the action-angle variables D^^ are an obvious choice. Note 
that in action- angle coordinates, the operator d/dcj)^ is invariant under the change 
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of origin while the operator d/dA^ depends on the origin (unless it is pulled 

out of an angle average, as in (|67|l '). Here we restrict attention to coordinates D^^, 
ignoring other possibilities, including various non-coordinate frames. The connection 
coefficients are required to vanish in coordinates so that the original definitions of 
diagrams remain unaltered (thus, the connection is fiat), but they take on non-zero 
values in the coordinates D'^. 

The standard (nontensorial) transformation law for the connection coefficients 
gives their values in coordinates D^, 

where in the second equality we have used 
^z^^ dz" 

- (70) 

to express one Jacobian in terms of its inverse. Only derivatives with respect to 
appear in the final expression (|69(l . The result is equivalent to 

dz d'^z" 

which is completely symmetric in indices (/Lta/3), as may be seen by differentiating 
®- . 

With these connection coefficients we can compute any diagram in terms of 
covariant derivatives in action- angle coordinates. For example, we find 

Ijo.^jjP^jj'y = Df^ J"" Df^^ J"^ = r"^T, (72) 

which as in (|fl.l|l shows the complete symmetry of this diagram in indices (a/37). 
Other diagrams can be expressed in terms of the connection coefficients, for example, 
by inserting resolutions of the identity we find 

^ £)/3 ^ ^ = r°^'" r^f"", (73) 

d'^^dp --r,"^r^''^ (74) 
/ \ =r^"T/^r,^^ (75) 

^ 

etc. As for the troublesome angle-dependent diagram in |(2Jl, it can be written in the 
interesting form, 

BY 

D'^^A^^D^ ^T^^--^LJ^, (76) 



where we use (|H.5|) . Note that F^""^ r^^^ = 0. The derivatives of the connection 
coefficients are not independent, due to the vanishing of the curvature tensor. 

4.4- Reduction to the one- dimensional case 

It is not obvious that the quantization condition ((SJ reduces to the one-dimcnsional 
formula (Fl) in the case = 1. To show that it does, we transform the 
multidimensional formula and then specialize to = 1. We begin with an identity, 

A. (A^- ^^k^A^^ A= ^ A^ ^A^^ 00 , (77) 
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where the second term on the right vanishes upon 0-averaging. We prove this by 
writing the left-hand side as 0' ^ {A^ ^ A'') , expanding this by the chain 

rule (instead of using l|H.5(l ^ to get three diagrams, of which two can be combined 
into the first term on the right hand side of (|77|l . transformed by a resolution of the 
identity. The third diagram of the expansion can be transformed by (|H.5|I into the 
second term on the right hand side of (|77|) . Now using Ht)8|) and (|77|l we can transform 
the order fi^ term of (0) (omitting the term iJ^ and the 1/48) into the form 

^A'^ A"') + 2u;jk,n [A'' ^ A' ^ A'") ) 




(78) 

The final diagram of this expression contains the dependence on the origin of the 
angles (before averaging). 

We now set = 1, dropping all Roman indices jkl, etc. Then the final diagram 
in (|78|l contains two identical (A) entries, and therefore vanishes. Then using the 
one-dimensional version of 177() . we can write the remainder of (|78|l as 

^{.iA^A)-2.'iA^A^A))^ = ^l-{H.^H,)^, (79) 

which agrees with the quantization rule cited in I if we expand, H = Hq + h?H2 + . . .. 



5. Conclusions 



The calculations of I and this paper achieved a normal form in a neighborhood of a 
stable fixed point of the Hamiltonian. Other authors (Helffer and Robert 1984, Colin 
de Verdiere 2004) have worked with an annular strip (in one-dimensional problems), 
that is, a range of energies bounded away from the ground state. This has the 
advantage that one does not need to worry about what happens inside the annular 
region: it may contain separatrices, nongeneric fixed points, etc., and the Bohr- 
Sommcrfcld rules are still valid. In fact, a WKB-Maslov kind of approach suggests 
that the asymptotic series for the eigenvalues can be computed in the immediate 
neighborhood of a single (quantized) torus. In this sense, our approach has been more 
global than it needs to be (for finding the eigenvalues). On the other hand, if one is 
studying normal forms, then a region containing a fixed point is certainly one of the 
cases to be studied. It is also a necessary case, if one wishes to divide phase space into 
regions, each of which can be transformed into some normal form. 

The normal form transformations we have considered (in one and many 
dimensions) are not the simplest. That would be the case of a contractible region 
of phase space in which the flow vector does not vanish, that is, one not containing a 
critical (fixed) point of the Hamiltonian. In this case the quantum normal form is the 
same as the classical, that is, the transformed Hamiltonian is just one of the momenta, 
say, pi. Integrability (a set of other operators that commute with H) is not required. 

As is well known, the torus quantization rule is not meaningful for nonintegrable 
systems, and neither is the normal form transformation of this paper. That is, in 
a typical (nonintegrable) system, in a neighborhood of a stable fixed point of the 
Hamiltonian, there is no classical normal form, and hence no quantum normal form 
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either, except in the sense of approximations vahd near the fixed point. On the other 
hand, the corresponding timc-dcpcndcnt problem, in which one uses an extended 
phase space with —H and t as conjugate coordinates, can be transformed into a 
simple normal form as described in the preceding paragraph. 

Conversely, instead of extending the phase space one can reduce it, in cases of 
symmetry. The case of SO{3) symmetry is especially important in practice (Littlejohn 
and Reinsch 1997). It would be interesting to study quantum normal forms on the 
reduced space, in which the wave function is a section of a bundle, and their relation 
to normal forms on the original phase space. 

Appendix A. Notation and terminology 

Wave functions in this paper are functions ip : M.^ C, that is, they represent 
a scalar particle in A^-dimensional space. The wavehmctions belong to the Hilbert 
space L^(R^), and operators such as act on this space. Operators are indicated 
with hats, and their Weyl transforms (or Weyl "symbols") are indicated by the same 
letter without the hat, for example, is the operator and its Weyl symbol. 

Latin indices i, j, etc. run from 1 to N. The superscript or subscript position of 
such indices, for example, , cuij, has no significance. Greek indices fj,, v, etc. run from 
1 to 2A^, for example, the coordinates on phase space are z'^ = (xi, . . . , x^^pi, . . . ,Pn)- 
Superscript and subscript positions of Greek indices are meaningful, and indices are 
raised and lowered with the Poisson tensor J'^" and its inverse J^w These are given 
in matrix form by 

"^"^ " ( -/ ) ' "^M- = ( ? ) ' (^•■^^ 

where / is the N x N identity matrix. Note the identity (I. A. 3). 

As in I, we use a nonstandard notation for mappings between sets. If A and B 
are sets, we write f : A ^ B to indicate that the domain of / is some suitably chosen 
subset of A (the standard notation makes A itself the domain). See Appendix I.B for 
notation concerning functions and their derivatives. 

In this paper, we denote the group of N x N integer matrices with integer inverse 
by GL{N, Z). The determinant of such a matrix is ±1. If the determinant is restricted 
to +1, we have the group SL{N,Z). 

Appendix B. Symbols and functions of harmonic oscillators 

In this appendix we outline the proof in the multidimensional case that an operator 
is a function of the set of harmonic oscillator operators { } if and only if its symbol 
is a function of the symbols {P} of these harmonic oscillators. 

The harmonic oscillators {/■'} commute and possess a simultaneous, nondegen- 
erate eigenbasis, the usual basis {|ni . . . rijv)} of harmonic oscillator states. Thus an 
operator (say, A) commutes with the operators {P} if and only if it is diagonal in 
this basis, which it does if and only if it is a fimction of the operators {/■'}. And 
A commutes with the {P} if and only if it commutes with the unitary operators 
lfi{t) =e^p{-itP/h), 

W{t)AW{ty =A, (B.I) 
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for all t. But according to Eq. (6.27) of Littlejohn (1986) (expressing the well known 
metaplectic covariance of the Weyl symbol, see Voros 1977, de Gosson 1997), when 
HB.1|I is mapped into Weyl symbols we have 

A{S^t)-^z) = A{z), (B.2) 

where S{t) is the classical linear symplectic map generated by P (that is, its 
Hamiltonian flow). But this map is equivalent to 9^ 9^ + t, showing that the 
symbol A{z) is constant on each torus of the set of classical harmonic oscillators {/-' }. 
This implies that the symbol A, regarded as a function of {9, 1), is, in fact, independent 
oi 0, A ~ A{I). Conversely, if A is such a function, then it satisfies HB.2|I and hence 
the operator A satisfies (jB.l|) and therefore commutes with the P . 



Appendix C. Normal forms for commuting, quadratic Hamiltonians 

The proof of the theorem (|48|l is a straightforward but somewhat lengthy exercise in 
linear algebra, so here we shall simply outline the highlights. The proof relies critically 
on the positive definiteness of the matrix (without this condition the conclusions 
do not hold). The notation of this appendix is independent of that of the rest of the 
paper, and will be explained as we proceed. 

Matrices and A-' are as defined in Sec. |21 In particular, the set {Q-'} of 
real, symmetric, 2N x 2N matrices is linearly independent, the matrices commute in 
the sense of ()45f) or (I46II . and is positive definite. Matrix J is the matrix with 
components J'^'^, defined in lA.ll and in this appendix we write uj = for its inverse 
(the symplectic form). Matrices and uj specify bilinear forms acting on 
C^^. We write 

for x,y £ or C^''^, where x is the complex conjugate of x, so that the symbols 
ui and are used both for the matrices and for the corresponding forms. Matrices 
JQ^, on the other hand, specify maps : M.^^ R^tv . (^2N _^ ^2N ^ 

The matrix JQ^ has purely imaginary eigenvalues that come in pairs ±jA, where 
A > 0. Let (Ai,...,Ar) be the set of distinct positive numbers that specify these 
eigenvalues. Let Ek C C^''^ be the eigenspace of JQ^ corresponding to eigenvalue iXk, 

Ek = {xeC^''\JQ^x^iXkx}, (C.2) 

and let Uk = dimEk- Then Ek is the eigenspace of JQ^ corresponding to eigenvalue 
~iXk- It is linearly independent of Ek, E^nEk = {0}, and dimii'fc = dimEk = Uk- Let 
= Ek®Ek C C^^, where the c means "complex," so that dimi^^ = 2nk- Subspace 
possesses a real basis, and so specifies a corresponding real subspace Fk C 
(the span of the real basis with real coefficients) , with dim Fk = 2nk ■ The subspaces 
Fk are symplectic vector spaces under the symplectic form oj, and thus possess a basis 
{gfci , ■ • ■ , qkuk , Pfci , ■ • ■ , Pknk } such that 

i^iqkhqkni) = i^{pkhPkm) = 0, 

^{qkUPkm) = -^^{Pkl,qkrn) = -Sim, 

for l,m = 1, . . . ,nk and k — 1, . . . ,r (Arnold 1989, p. 219). Moreover, 

Fi®...®Fr= R^^, (C.4) 



(C.6) 



(C.7) 
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so that 

r 

Y^nu^N. (C.5) 

k=l 

This means that the matrix JQ^ can be diagonaUzed by a Hnear transformation, that 
is, its Jordan normal form has no off-diagonal elements. 

The {<;,p} bases on the subspaces Fk can be chosen so that the matrix elements 
of are a multiple of the identity on each subspace, 

Q^iqkUqkm) = Q^iPkhPkm) = XkShn, 
Q^{qkUPkrn) = Q^{PkUqkm) = 0, 

and so that the matrix elements of the other are diagonal in the form, 

Q^iqkhqkm) = Q^iPkhPkrn) = al^Slm, 

Q^qkhPkm) = Q^ipkhqkm) = 0, 

for j > 1 (no sum on I), where the a^^ are real numbers. 

We create a single index for the q and p basis vectors by writing 

qt = {qi, ■ • • ,gjv) = {qn, ■ ■ ■ . --qri, ■ ■ • 

Pi = {Pi, ■ • ■ ,Pn) = (Pii, • ■ • ■ • -Prl, • ■ • ,P™J, 

for i — 1, . . . , TV, and similarly we write 

= (-^1: ■ • ■ I '^l, ■ • ■ 7 ^r, ■ ■ ■ , Ar), (C.9) 

with each Afc repeated times, and 

a^i = (a{i,...,aj„^,...,a^i,...,a^„J, (C.IO) 

for J > 1. We also write an = A^. We arrange the basis vectors columnwise in the 
order {q,p) to create a 2N x 2N matrix S. Then ljC.3|l is equivalent to S'^ujS — lo, 
so that S is symplectic, and (|C.6p and HC.7|I are equivalent to which was to be 
proven. The N x N matrix is nonsingular because the matrices are linearly 
independent. 



Appendix D. An example of the classical normal form transformation 

In this appendix we consider central force motion in two dimensions, in which the 
Hamiltonian is 

H{^,p)^^ + V{r), (D.l) 

where x = (x, y), p = {px,Py), r = |x| and p — |p|. We omit the hats, but H can be 
either the quantum Hamiltonian or its Weyl symbol (the classical Hamiltonian). The 
"other Hamiltonian" is the z-component of the angular momentum, 

L(x, p) = xpy - yp^. (D.2) 

The Weyl symbols of these two operators have no h dependence, so the principal 
symbol and the whole symbol are identical, and we shall simply write Hj (with a 
subscript) instead of Hq. Thus, Hi = H and H2 = L. 

We assume the potential is smooth and has a minimum at r = 0, so that x = p = 
is a stable fixed point of the classical motion. We assume V^(0) = so that H and 
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L vanish at the fixed point. We write Wo = y^V"{0)/m for the frequency of smaU 

1 /2 

vibrations near the origin. We transform to variables with balanced units of action ' 
by writing x = x'/y^mu^, p = y'mujo p', substituting these into the expressions for 
H and L, and then dropping the primes. The function L does not change under this 
transformation, but H becomes 



2 



2U{r)], 



(D.3) 



where U{r) {1 / ujo)V {r / ^/rnuJ^) . Note that U"{Q) = 1. We transform (jn3|l to polar 
coordinates {r,9), to obtain 

L2 



H = 



UJQ 

2 



Pr 



2U{r) 



(D.4) 



where Pr = x • p/r and L — pg. For positive energies E below the first separatrix (if 
there is one), the classical motion in the x-y plane takes place in an annular region 
between two turning points ro, ri. 

This region is the projection of the 2-torus in the 4-dimensional phase space onto 
the x-y plane. The torus is characterized by the values of the energy E and angular 
momentum L (in the latter case we confuse a function with a value of a function). 
The turning points are the roots of (jD.4|) obtained by setting pr — 0, that is, the roots 
of 



^2 

2E/UJ0 = — 



2U{r). 



(D.5) 



There are two branches to the projection of the torus, identified by the sign of pr. The 
branches are glued together at the radial turning points, where Pr = 0. The angular 
momentum L can take on positive or negative values. The inner turning point ro goes 
to zero as L passes through zero for fixed E, so the projection of the torus closes up 
in the center and then opens up again. Exactly at L = 0, the projection of the torus 
onto the x-y plane is singular at the origin (the inverse image of the origin is a circle 
on the torus in phase space). 

To study the normal form transformation for this system in a small neighborhood 
of the fixed point, we expand (|D.3(I for small r to obtain, 

^0 / 2 , 2\ ^0 , 



y 



Pv 



(D.6) 



From this and l|D.2p the matrices and defined by (|41|l are easily read off. In 
particular, = uiqI (/ is the 4x4 identity matrix), so is positive definite. 

The symplectic matrix S in H48|l can be found by various means, for example 
by noting that in this problem quadratic Hamiltonians H and L belong to the Lie 
algebra of the U{2) subgroup of S'p(4), and using the theory of SU{2) rotations to 
transform L into a diagonal form. The U{N) subgroup of Sp{2N) is basic to the 
theory of the isotropic iV-dimensional harmonic oscillator (Littlejohn 1986, Cushman 
and Bates 1997). We omit the details and just quote the result, which we write as a 
linear canonical transformation. 



(D.7) 



/ X \ 












-1 








y 


1 





1 


-1 











Px 


V2 





1 


1 









Pi 


\ Pv 1 




1 1 








1 


) 




\ V2 ) 



Multidimensional quantum normal forms, Moyal star product, and torus quantization26 



where the matrix shown is the matrix S of H48|l . When HD.6|I and HD.2|I are transformed 
to the new coordinates, we find 



(D. 



HoZ-' = -{xl+pl + xl+pl)=uj„{h+l2), 

LoZ-' = \{xl+pl-xl-pl)^h-h, 
where /i, I2 are defined by ((SJ- From this we can read off the matrix Ojk of l(5 

The actions Aj depend on the choice of basis contours on the tori. The most 
obvious choice is (rriFg), the contours whose projections onto the x-y plane take 
place at constant 9 and r, respectively. These are the contours that emerge from the 
standard method of separating the Hamilton- Jacobi equation in polar coordinates. 
Contour goes the "short way" around the torus, that is, from the inner to the outer 
radial turning point and back again, with pr positive on the way out and negative 
on the return. Contour Tg projects onto a circle in the x-y plane, traversed in a 
counterclockwise direction, with fixed value of {r,pr). The action associated with Tq 
is the angular momentum itself, 

Ag^^jpgde^L, (D.IO) 

while the radial action is 

Ar^ — iprdr^- J ^-2U{r)dr. (D.ll) 

2n J TT y LUo r^ 

A tricky aspect about these contours and actions is that the contour F^ is not 
continuous as L passes through zero for fixed E, that is, F,. for small positive L is not 
close (in phase space) to the contour F^ for small negative L. Instead, the contour 
at small negative L that is close to F^ for small positive L is F^ — Fg, that is, the 
contour, continuously tracked as L passes through zero from positive values, picks up 
an extra (backwards) loop around the origin. As a result, (F^, Tg) and (A^, Ag) do not 
constitute a smooth assignment of basis contours and actions in the region of interest, 
but rather only over the two subregions i > and i < separately. 

In the quadratic approximation (|D.6|) to the Hamiltonian (an isotropic harmonic 
oscillator) the radial action Ar can be evaluated explicitly. In this case the energy and 
angular momentum obey the bounds —E/luq < L < +E/uJo, and the radial turning 
points are 



ro,i = ^Je/luo ± ^{E/uj^Y-L-^. (D.12) 
Then direct integration of (|D.11|) gives 

Ar = \(--\L\). (D.13) 

The absolute value signs in the second term is an indication that F^ and hence Ar are 
not continuous at i = 0. 

For the quadratic approximations, a different contour on the tori is conceptually 
useful, namely, the one generated by H itself. This is possible because the isotropic 
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harmonic oscillator possesses only periodic orbits. Call this contour T^. By direct 
examination of the elliptical orbit in the x-y plane, it is obvious that 

Th = 2rr±rg, (D.14) 

where the ± is the sign of L. Thus, 

Ah ^ (f p ■ dx = 2Ar + \L\. (D.15) 



2vr jr. 



The contour Th is continuous as L passes through zero; one function of the absolute 



value sign in ljD.15|l is to cancel the discontinuity in the term 2Ar. The action Ah is 
easily evaluated in rectangular coordinates with time as the parameter of the loop by 
using Hamilton's equations of motion. The result is 

Ah^— p-xdt = E/LJn. (D.16) 

271- Jo 

This integral is much easier than the one that led to ljD.13|l . but combined with (|D.15p 
it gives the same result for Ar. 

The orbits of the Hamiltonian H are closed curves only in the quadratic 
approximation (|D.6|I . not for the full Hamiltonian (|D.3|) . Nevertheless, the closed 
contour Th, defined on tori near the fixed point, can be extended continously to 
closed contours on tori at larger values of E, where (jD.14|l and hence (|D.15|I are still 
valid. 

Let us choose the basis contours to be (ri,r2) = (ri.,r6)), so that {Ai,A2) = 
{Ar, L). Then contours (Ci, C2), the basis contours on the tori of harmonic oscillators 
(Ji, J2), when mapped by Z back onto contours on the tori of {H,L), can be worked 
out in the {x,y,Px,Py) coordinates by using (|D.7|) . For example, to get Ci, we may 
write xi — \J21x sin 6*1, p\ = \J2I\ cos^i, xi = 0, p2 = V2I2, substitute into (|D.7|) . 
and let Q\ range from to 27r. By plotting these contours in the x-y plane and keeping 
track of the sign oipr, we obtain the following relations among contours: 



i>0; '■' ^ L<0, (D.17) 





where T^ and Tq really mean Z oT^ and Z oTq. Integrating the symplectic form over 
these contours, we obtain. 



L>0; /' ^ L<0, (D.18) 





where 1\ and li really mean l\o Z and I2 o Z. These relations may also be obtained 
(more easily) by combining (|D.13|) and ljD.8|) . From them it is easy to obtain the 
matrix i/jk of (|55|) . We find, 

I \ y L>0; ,.^(^1 L<0. (D.19) 

Appendix E. Existence of the generators G2n 

The notation of this appendix differs slightly from that of the main text, in that we 
write for some intermediate step in the second stage transformations, 

= t/2„-2 . . . UiU2K'uM ■ ■ ■ uL-2, (E.l) 
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where n > 1 (rather than the final product as n — > oo, as in the main text). We also 
write 

M'^ = U2n M' UL (E.2) 

for the next stage in the transformation. Symbols and M'^ of and M'^ are 
expanded in H as in 1^5)1 . We assume that Mq = and that M2„, expressed as a 
function of 9 and /, is independent of 6 for < to < n — 1. Under these assumptions 
we will show that exists such that = for < m < n — 1 and such that 
Mj'Jj is independent of 9 (that is, to two higher orders in h than the previous step). 
Since the stated assumptions hold for n = 1, by induction they hold for any n > 1. 

The unitary operator t/2n is given in terms of its generator G2n by (1.47), where 
here we assume that G2n is independent of h. Thus, (|E.2(I has the expansion. 



M'^ = W - in^^-^[G2n,M'] - lh^"-^[G2n, [G2n, M']] + ■■■■ (E.3) 

Upon transcribing to symbols and expanding as in H25() we find Afjm — -^^2jn ^'^^ 
< m < n ^ I, and 

Ml = Ml + {G2n, M^} = Ml + ^,k^. (E.4) 

Taking the averaged part of this wc obtain Ml = Ml, that is. Ml is independent 
of 9. Taking the oscillatory part we obtain a differential equation for G2m 

^-^ = -{^-%.Ml. (E.5) 
A solution for this exists if the second derivatives 

are symmetric in fc and I, where we drop the tilde on Ml since the 9^ derivative of 



the averaged part vanishes. Rearranging, we see that a solution G2n of (|E.5|I exists if 
and only if 

dMl dMl ^ 

-c.,^+..,^=0, (E.7) 



that is, if 



{M^,Ml} + {Ml,M^}^0. (E.8) 

However, since the operators M^ are unitarily equivalent to the and since 
the latter commute, we have [M'^,M'] = 0. Transcribing this to symbols and using 
(I. A. 6), we obtain the hierarchy of equations, 

{MlM^} = 0, 

{Ml Mi) + {Ml M^} - l{Mo^ M^}3 - 0, 

{Ml Mi) + {Ml Ml + {M| , A/^} ~ ^{Mo^ M^jg 

-Y^iMlMl,} + J^{Mo^ M^}5 = 0, (E.9) 

etc. 

Now let A and B be any two functions of / that are independent of h. Then 
all the odd order Moyal brackets {A, B^2n+i vanish. This follows from the fact that 
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the corresponding operators A and B are functions of / (see [Appendix Bt , and hence 
commute, [A, B] = 0. Expanding the symbol version of this commutator as in (I. A. 6) 
then yields {A, B}2n+i = for all n > 0. 

Applying this theorem to (|E.9|) and using the fact that Af|^ is independent of 9 
for 0<m<n— l,we see that all terms in the first n — 1 equations of the hierarchy 
vanish, while the n-th one gives precisely the integrability condition (|E.8(I . Thus, G2n 
exists. 



Appendix F. An antiderivative for ()39|1 

We wish to evaluate the integral 

f de (G.^Ai) oZ-\ (F.l) 
Jo 

to within terms that vanish on taking a 0-average (not a (/)-average, as in %V!\ , because 
we have removed the oZ seen in H39|l 'l. We will write the integrand in the form 
{d/de){F^ o Z~^) for some function to be determined. According to (1.24) and 
(1.25), satisfies. 

This is the sense in which we shall regard F^ as the antiderivative of the integrand of 

For notational convenience we now drop the e-subscripts on A^, 0^, Z?^, Z^, F^ 
and Ge until further notice. That is, A^ , (fP , D^, Z, etc., now stand for quantities 
evaluated at any e (not just e = 1, as in the main text or in Table 

The functions F we are interested in are made up of diagrams involving the 
action-angle variables D^^. The derivative dF/de is evaluated by the chain rule, in 
which 

^ ^ i^E^oZ) = (Ell o Z)— = {Ei^ o Z){Z^ G} = {E^ o Z, G} 
de de de 

= {D^',G} = D^^G , (F.3) 

where we use H2()|l in the first equality, the chain rule in the second, (1.19) in the third, 
and the chain rule again in the fourth. See the first three terms of (jF.7(l for an example 
of the application of this formula. 

Thus, taking the derivative D /De of a diagram creates a series of new diagrams, 
each with one extra arrow and one extra copy of G compared to the original diagram. 
Since the integrand of Hi^'.l() has one G, one copy of A^ and three arrows, the suggestion 
is that the antiderivative must be a diagram with A^ and two arrows. Unfortunately, 
this cannot be correct, since any diagram with an arrow that begins and terminates on 
the same operand {A^ in this case) vanishes, by the antisymmetry of J^". The basic 
logic is correct, however, it is just that we must take into account the possibility of 
introducing a resolution of the identity (see [Appendix C\ into the original integrand, 
which adds one copy each of D^^ and I?^, and one extra arrow. Thus, the original 
integrand can be regarded as a diagram with I?'', Z?^, A^ , G, and four arrows. 
Therefore candidates for the antiderivative would involve diagrams with D^^, D^, A^ 
and three arrows (and no G). It would also be possible to introduce more than one 
resolution of the identity, but it turns out that one suffices for our purposes. 
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Therefore candidate diagrams for the antiderivative include the foUowing: 

AJ^D^^D^, / \ , (F.4) 

£)P — ^Df, 

where we do not count as distinct any diagrams that can be obtained from the ones 
given by reversing the direction of arrows or by using (I. A. 3) . Let us call these diagrams 
Ti, T2, T3 and T4, in the order given. 

These diagrams are not all independent. First we note that the Poisson bracket 
{D^^, D'^} = J^'^ is a constant, so any derivatives of it vanish. The action A-' can 
be regarded as a component of some 2 A^- vector , so the same applies to Poisson 
brackets {D^',A'}, etc. 

These constant Poisson brackets lead to various identities, such as 

= (^^ ^ ^m) = ^ -^D^ +2 / \ 

+ Df'^D^^ = 2(T3 + T4), (F.5) 



and 



0= Df,^ {A^ D^) ^ D^^A^ ^ + 2 / \ 

A3 D^' 

+ A' ^ D'' 5= = -Ti + 2r4 + T2, (F.6) 

where we expand using the chain rule for diagrams. We see that only two of the 
diagrams HF.4|) are independent, say, Ti and T2. 

We take the first of these Ti as a guess for the desired antiderivative, and compute, 

R. (i^M ^ ^ D^) = [D^^ ^ G) ^ A^ z^D^ + [A^ -^G)^D^ 

+ A^ ^ {Df' ^G) - {Df" A^ i^D^)^G , (F.7) 

usmg (|R2)| and (jESl- The four dia grams on the right hand side can be expanded by 
the chain rule, resulting in fourteen diagrams altogether, but ten of these cancel just 
by changing the directions of arrows or using (I. A. 3). The remaining four diagrams 
add up to a continuation of the previous equation, 

t X tt 
= jjt.^ A^ ^G^D^ - 2D^^G^A^ _ D^'^G^A^ 

_ ^ =>: G ^ D^, . (F.8) 

We notice the following facts about this result. First, the final diagram is the 
one we want, G ^ A^ ■, slightly disguised with a resolution of the identity. Next, we 
are allowed to split off an exact (/>-derivative (i.e., a -derivative, which turns into a 
^-derivative when composed with Z'^^). But if X is any diagram, then 

rlX 

A^^iX)=-^, (F.9) 
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as we see by expanding the Poisson bracket in ((/>, A) coordinates. So to create an exact 
^-^ -derivative, we must create diagrams of the form {X) . Third, to create such 

diagrams we can employ identities hke ljF.51) and IjF.Gp . which follow by using some 
constant Poisson bracket as a part of a larger diagram. In the present case, useful 
identities follow from the vanishing expressions. 



G: 



G: 



(D^^A^). (F.IO) 



{D^' 



G^D^),{FA1) 



By using these we obtain finally, 

§-jD^^A^^D,)=2G^A^-^^ 

showing that the first guess for the desired antiderivative worked (the other one, the 
second term in (jF.4|) . called T2 above, does not yield the desired integrand). 
Now restoring the e's and using (|F.2|I . we have 



de (Ge 



de 



Al) oZ-'^ (D^ 
d 



■Al 



:D„) oZ- 



oz: 



P 



E„ 



a 



(F.12) 



The middle term on the right hand side vanishes, since P is a quadratic polynomial 
in z'^ with three arrows attached to it, and in the final integral we have used the 
composition with Z^^ to convert the 0^ derivative into a 9^ derivative, which can be 
taken out of the integral. Dividing the result by 2 and composing with Z gives (|39|) . 



Appendix G. Resolution of the identity 

The technique we call "inserting a resolution of the identity" can be used to break 
a bond of a diagram, resulting in a new diagram (or product of diagrams) with one 
extra arrow and an extra copy of Z?^ and D^. In this appendix. A, B and C are any 
functions, while D^^ and are as in the main text (see (@J). 

Consider a diagram (or product of diagrams) of the form X , where 

X stands for a single function, a diagram with more than one function connected 
by arrows (where the arrows coming into or out of X are attached anywhere on 
the diagram), or a product of such diagrams. Examples are i'/i — ^ A — *- Z)^ , 



A^B^D^' , or {D. 



D,. 



X^Di' 



■A) {B 
dD 



(7 ^ , Then we have 



Jl^" „Xr J^" 



(G.l) 



But since the components of the symplectic form are preserved under a canonical 
coordinate transformation, we have 



dz°' dz^ dz 



dzP 



J, 



Q/3, 



(G.2) 



Multidimensional quantum normal forms, Moyal star product, and torus quantization32 
and therefore, 



D^^X^D^ =r'' .Xr^ ^X-^ _ (G.3) 



For example, 



D^^A^Dt' - ^A-^ =0, (G.4) 

D^,^ Df" ^ A^B =^ - A^B , (G.5) 

(D^^A) {B^C^D^') = B^C^A. (G.6) 

Appendix H. Diagrammatic identities 

In this appendix we present some diagrammatic identities involving the action-angle 
variables D'^. 

For a linear diagram with three D's, we have 

D^^ ^ D" D"" = symmetric in (^j/a) . (H.l) 

This follows by expanding vanishing diagrams such as {D^^ D^) D" by the 
chain rule, and by inverting the directions of all arrows. 
Similarly, the triangle diagram, 

(H.2) 

/ \ 
L>° ^ 

can be shown to be symmetric (antisymmetric) under even (odd) permutations of 
(a/37). 

For a linear diagram with four £)'s, let us abbreviate 

D^^D'^^D^^D^' = (H.3) 
Then we have the symmetries, 

[^va(3] = [v^aP] = [^vl3a\ = — [aP^v], (H.4) 



which follow by the same techniques used in proving (|H.ip . 
Finally, we have the identity, 

£)M ^ (D" ^ ^ 1)7) ^ ^ . (H.5) 

We prove this by directly expanding the diagram on the left by the chain rule and 
then using the expansion of the vanishing diagram ^ [D^^ D^) and 
the identities ljH.411 . Note that both sides of IjH.Sp are manifestly symmetric in (a/37) • 



Multidimensional quantum normal forms, Moyal star product, and torus quantizations^ 
References 

Abraham R and Marsden J E 1978 Foundations of Mechanics (Reading, Massachusetts; 

Benj amin / Cummings) 
Alber M S, Luther G G, Marsden J E, and Robbing J M 1998 Physica 123D 271 
Arnold V I 1989 Mathematical Methods of Classical Mechanics (New York; Springer- Verlag) 
Bayen F, Flato M, Fronsdal C, Lichnerowicz A and Sternheimer D 1978 Ann, Phys, 111 61 
Berry M V and Tabor M 1976 Proc. Roy. Soc. Land. A 349 101 
Brillouin M L 1926 J. Phys. 7 353 

Cargo Matthew, Gracia-Saz Alfonso, Littlojohn R G, Reinsch M W and Rios P de M 2004 J. Phys. 
A 38 1977 

Charbonnel, A-M 1986 Integral Equations Oper. Theory 9 502 
Child M S 1998 J. Phys. A 31 657 

Colin de Verdiere Y 2004 Bohr-Sommerfe ld rules to all orders Preprint 

http: / / www-fourier . uj f-grenoble. fr~ycolver / ebk. ps' 
Cushman R H and Bates L 1997 Global Aspects of Classical Integrable Systems (Basel; Birkhauser 

Verlag) 

Cushman R H and Duistermaat J J 1988 Bull. Am. Math. Soc. 19 475 

Cushman R H, DuUin H R, Giacobbe A, Holm D D, Joyeux M, Lynch P, Sadovskii D A, and Zhilinskif 

B I 2004 Phys. Rev. Lett. 93 024302-1 
de Gosson M 1997 Maslov Classes, Metaplectic Representation and Lagrangian Quantization (New 

York; Wiley) 

Duistermaat J J 1980 Commun. Pure Appl. Math. 33 687 
Eckhardt B and Wintgen D J. Phys. A 24 4335 
Einstein A 1917 Verh. dt. Phys. Ges. 19 82 
Eliasson L H 1990 Comment. Math. Helvetici 65 4 
Fedosov B V 1994 J. Differential Geometry 40 213 

Fomenko A T 1991 Advances in Soviet Mathematics vol 6 (Providence, Rhode Island; American 

Mathematical Society) 
Foxman J A and Robbins J M 2004 preprint m ath-ph/0411017, 
Gracia-Saz A 2005 preprint math.QA/0411163 

Helffer B and Robert D 1984 Ann. Inst. H. Poincare Phys. Theor. 41 291 
Keller J B 1958 Ann. Phys. 4 180 

Kirillov A A 1976 Elements of the Theory of Representations (New York; Springer- Verlag) 

Littlejohn R G 1986 Phys. Reports 138 193 

Littlejohn R G and Reinsch M 1997 Rev. Mod. Phys. 69 213 

Littlejohn R G and Robbins J M 1987 Phys. Rev. A 36 2953 

Marsden J E and Ratiu T S 1999 Introduction to Mechanics and Symmetry 2nd ed (New York; 
Springer- Verlag) 

Maslov V P and Fedoriuk M V 1981 Semi- Classical Approximations in Quantum Mechanics 

(Dordrecht; D. Reidel) 
Percival I C 1973 J. Phys. B 6 L229 
Voros A 1977 Ann. Inst. Henri Poincare 4 343 

Weinstein A 1977 Lectures on Symplectic Manifolds, Regional Conference Series in Mathematics 29 
(Providence, Rhode Island; American Mathematical Society). 



Multidimensional quantum normal forms, Moyal star product, and torus quantizationZA 
Figure captions 



Figure 1. The two energy-momentum maps (the -ffo-map and the /-map) are 
related by map Z and the matrix a. The figure illustrates the case N = 2. 



Figure 2. Maps associated with the region of interest, R. The four axes in the 
diagrams on the left represent IR^'^ , and the two axes in the diagrams on the right 
represent IR''. 



